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ABSTRACT 

We present a coanand generator tracker approach to model following control 
of linear distributed parameter systems (BPS) whose dynamics are described on 
infinite-dimensional Hilbert spaces. This method generates finite-dimensional 
controllers capable of exponentially staple tracking of the reference trajector- 
ies when certain ideal trajectories are known -to exist for the open-loop DPS; 
we present conditions for the existence of these ideal trajectories. An adaptive 
version of this type of controller is also presented and shown to achieve (in 
some cases, asymptotically) stable finite-dimensional control of the infinite- 
dimensional DPS. 


I. INTRODUCTION 

By a distributed parameter system (DPS) , we mean a system whose dynamical 
behavior with respect to external disturbances is described by partial differ- 
ential equations. Of course, everything is a DPS if it is carefully scrutinized, 
especially if high performance is demanded, e.g., a simple electrical circuit 
at very high frequencies. However, lumped parameter (ordinary differential 
equation) approximations often suffice to describe the system behavior of many 
engineering systems. Indeed, such approximations are necessary for DPS control- 
ler designs to be Implemented with on-line digital computers. Nevertheless, 
the distributed parameter nature of control problems should not be discarded 
prematurely; otherwise, control approaches can be generated which look good on 
paper but are not sufficiently robust to operate with the actual system. This 
has been illustrated in computer simulation and in even a few laboratory demon- 
strations of flexible structures, yet, it continues to be ignored in some parts 
of the control community. To understand the controller-structure interaction, 
a DPS viewpoint is essential . 

The most rerious difficulty of the DPS viewpoint is that it requires the 
mathematical ideas of infinite-dimensional function spaces and unbounded oper- 
ators on these spaces; for example, see [1]— [2]. Several results in the past 
have been posed within this mathematical framework with the required mathematical 
rigor [3]. Yet, the necessary practical constraints were interpreted so that 
the results would be relevant to structural dynamiclsts and control system 
engineers and would make the maximum use of their experience and intuition. 

With these ideas in mind, the concept of model following appears to be a 
procedure that yields a useful finite dimensional controller that might be 
designed taking into account the distributed nature of the system dynamics, 
whereas early model following control systems required the satisfaction of 
certain ,: Perfeet Model Following" conditions which necessitated the use of a 
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reference model having the same order as that of the process [A] , the more 
recent output model following controller or Command Generator Tracker (CGT) as 
developed by Broussard [51 allows the use of a model of arbitrary order, provided 
that the number of controls is equal to the number of outputs being controlled. 
This concept in fact served as the basis for a finite dimensional adaptive 
controller that was used for controlling large structural systems [6, 7]. 

Thus since the CGT algorithm makes it possible to use a finite dimensional 
reference model which subsequently gives a finite dimensional controller regard- 
less of the process order. This provides the basis for a direct adaptive con- 
troller which produces stable closed-loop operation with the class of linear 
distributed parameter systems considered here. The difficulties of stable 
adaptive distributed parameter control are detailed in, e.g., [8]-[9] and the 
references contained therein. In Sections 2 and 3 the nonadaptive model 
following controller is developed and analyzed; in Section A, the adaptive 
version is presented and shown to produce a stable closed-loop. Conclusions 
and future directions are presented in Section 5. 

2. PROBLEM FORMULATION 


2.1 Process Description 


The distributed parameter systems (DPS) of interest will be modeled by the 
following state space form : 


j “ Av(t) + Bf (t) ; v(o) = v o 

1 y(t) = Cv(t) 


(2.1a) 

(2.1b) 


where the state v(t) is in an infinite-dimensional 
inner product (•, •) and corresponding norm 1 1 • j | . 
operators B and C have the sam e finite rank P, and 
i nputs for P linear actuators and the outputs from 
ly. Thus, 

P 

Bf (t) = l b f.(t) 

1=1 


real Hilbert space H with 
The bounded input-output 
f(t), y(t) represent the 
P linear sensors, respective- 


( 2 . 2 ) 


and 


y(t) = [y 1 (c). 


v p (t)] 


with 


y^(t) = (Cj , v(t)) ; 1 <_ j < P 


(2.3) 


where and belong to H. 


In infinite-dimensional theory, the operator A 


is a closed, linear, unbounded (differential) operator with domain D(A) dense in 
H. Furthermore, (2.1)-(2.3) represents some well-posed physical system, which 
in mathematical terms is the weak formulation of (2.1): 
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(2.4) 


t 

v (t ) = U(t) v q + / U(t-r) Bf(T)dT 
y(t) = Cv(t) i t > 0 

where is any initial state in H and U(t) is the C^-semigroup of bounded 
operators generated on H by A, This latter means: 


U(t+x) = U(t) U(x) ; 

t > 0, T >_ 0 

(2.5a) 

U(o) = I 


12.5b) 

lim (U(t) - I] v = 0 
t-K) + 

; v in H 

(2.5c) 

Av = [lim ^ )v ; 

v in D(A) 

(2.5d) 


t-0 + 


Note that the semigroup U(t) evolves the initial conditions v^ forward in time. 

When v is in D(A) and f(t) has continuous first derivative, v(t) also is differ- 
entiab?e, lies in D(A) for t >_ 0, and satisfies (° 1). However, any v q and H 

and any square-integrable f(t) wili satisfy the weak formulation (2.4) and yield 
states v(t) in H for all t ^ 0. Consequently, (2.4) is easier to work with in 
infinite-dimensions and is more likely to represent the actual physical system 
being modeled by (2.1). This form, (2.1) or (2.4), models most practical 
interior control problems for linear DPS where the actuator and sensor influence 
functions are given by b^ and c , respectively. 


For example, control of the damped wave equation on a region OCR 
single actuator and sensor is described by (for e > 0) : 


by a 


t-aiflL ♦ t . auu.o . b(l) f(t ) 

dx 

y (t) = / c(x) u(x,t) dx 

a 


(2.6a) 


(2.6b) 


where u(x,t) is the displacement from equilibrium of ft and the influence func- 
tions b and c can be taken as approximations of Dirac delta functions at the 
location of the actuator and sensor. 


A u(x,t) = l 
° l«l 


3 u(x,t) 
3x 2 „ 


The operator A q is the Laplacian given by 

(2.7) 


on D(A ) = {u(x,t)eH ju(x,t) is smooth and u(x,t) ' 0 on the boundary of Q}. 

O O A 

The domain D(A ) is dense in H = Ir(fl) with the usual inner product (♦,*). 

o o ° 

This can be put into the form (2.1) by choosing the state v(t) * {u(x,t). 
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SuU.Q .T 

3t J 


in H = D(A 1/2 ) 
o 


x with the energy inner product : 


(v, U ) = (A c 1/2 v r A o l/2 Wl ) o + (v 2 , « 2 ) o 


( 2 . 8 ) 


The operator A in (2.1) becomes 


A = 



I 

-el 


and the rest follows. 


(2.9) 


Another important example is the mathematical setting for large structural 
systems (LSS) which may be described as a continuin by the following system of 
partial differential equations: 


m(x)u tt (x,t) + D Q U t (x,t) + A Q u(x,t) - F(x,t) (2.10) 

where u(x,t) represents a vector of instantaneous displacements of the structure 
0 from its equilibrium position due to transient disturbances and the applied 
force distribution F(x,t). The displacements can be translational and rotation- 
al, and the forces can be generalized to include torques, as well. The mass 
density m(x) is positive and bounded on ft. 


The internal restoring force term A q u is generated by a time-invariant, 
symmetric, non-negative differential operator A^ appropriate to the LSS. The 

domain D(A ) of A contains all smooth functions satisfying the LSS boundary 
conditions°and is°dense in the infinite-dimensional Hilbert space * L (ft) 

with the usual inner product (*,*) o and associated norm ||*|| o . In most cases, 

the operator is assumed to have discrete spectrum, i.e.. Isolated resonances; 

this can be expressed by the following eigen-problem: 

Vk = w k *k (2 - n) 

where w^ are the vibration mode frequencies and are the corresponding 

vibration mode shapes . Of course, exact expressions for this modal data are 
rarely known for an actual LSS . 

The damping term D Q u t is composed of a skew syvnetrlc part, which repre- 
sents gyroscopic damping due to any on-board rotors or constant spin rate of 
the whole LSS, and a small symmetric part which represents the internal struc- 
tural damping and is thought to provide very low mode damping. 

The applied force distribution is 

F(x, t) * F (x, t) + F (x, t) (2.12) 

c D 
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where F Q represents the external disturbance forces on the LSS (and possible 
nonlinearities) and represents the control forces due to P actuators: 


F 

c 


B f 
o 


P 


l 

i-1 


b i (x) f i (t) 


(2.13) 


where the actuator amplitudes are f^t) and the actuator influence functions 
are b^x) in H q . These are usually localized or point devices so that they 
approximate fi(x-x^); however, they do not have to be point devices. 

Observations are obtained by P sensors 


y = c 0 u + E 0 u t (2.14) 

where y.(t) = (c^ ,u^) + (ej,u t ) o , 1 <_ j P, with influence functions c^ for 
position sensors and e^ for velocity sensors in H q . Again, these are usually 
localized or point devices but they do not have co be. 

The LSS dynamics are defined by (2.10) and (2.14) can be put into the 
inf inite-dinensional state space form : 

| - Av(t) + Bf (t) + rf D (t) (2.15a) 

i y(t) = Cv(t) ; v (o) = v o (2.15b) 

with (A,B,C) as in (2.3) and the persistent disturbance term Tf p (t) obtained 

from F Q in (2.12). Impulsive disturbances in the structure are modeled by the 

initial condition v . 

o 

The Kille-Yosida Theorem (e.g. [1], Theo. 8, 9, p. 153), provides condi- 
tions under which an operator A generates a (^-semigroup U(t) satisfying: 

| |u(t) | | £ Ke _0t , t > 0 (2.16) 


where K >_ 1 and o real. The necessary ^nd sufficient conditions are given for 
the resolvent operator R(X,A) = (XI-A) : 

| ! R(X , A) n | | < - K — ; n = 1, 2, ... (2.17) 

(X+o) n 

for all real X > - o in the resolvent set of A, p (A) * £X complex |r(X,A) is a 
bounded operator on H} . The spectrum of A, o(A) = p(A) is much more compli- 
cated in infinite-dimensions, but, in finite-dimensions, it consists only of 
the (finite number of) eigenvalues of A. We say that A is exponentially stable 
when o > 0 in (2.16), i.e., the semigroup U(t) generated by A decays exponen- 
tially at the rate a. There are many other types of stability in infinite- 
dimensions, but no others provide the safety of a stability margi.i o; therefore. 
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this seems to be the kind of stability of most practical interest for engineer- 
ing applications where there is always some uncertainty in the model of JPS. 

2.2 Model Following Control Problem Formulation 

Given the DPS as defined in (2.1), it is desired to find a finite dimen- 
sional controller so that the output y(t) "follows" a desirable output trajec- 
tory y m ( t - This output trajectory is to be generated by the finite dimensional 

(asymptotically) stable reference model: 

« = V + Vm (2 ‘ 18a) 

= c B q ; ^<°) = % (2.i8o) 

where 

q is the model state vector having dimension N, 

u q is a step or reference level command with dimension P, 

y^ is the output trajectory also having the dimension P, 

and A , B are matrices with appropriate dimensions. It should be noted that 
mm 

the dimension of both y and u is the same as the dimension of the process 

to m 

input f and the process output y as defined in (2.1). Usually q * 0 will be 
chosen. 

The output model following control problem to be solved is the development 
of an algorithm that defines the process input f(t) so that the following two 
model following conditions (MFC) are satisfied: 

MFC 1) If yCtp = y m (t 1 ) , then 

y(t) = y m (t) . f »r t > t 1 

MFC 2) If y(t.) 4 y (t.), then 

y(t) asymptotically will approach y m (t), i.e. 

Aim [y(t) - y (t)] - 0 

t-x» 

3. DEVELOPMENT OF THE NONADAPTIVE MODEL FOLLOWING CONTROLLER 
3.1 Solution Definition 

In a manner similar to Broussard's development of the Command Generator 
Tracker (CGT^ [5], the concept of an ideal state v , control f and output 
trajectory y will be introduced. It is required that these trajectories 
satisfy the process dynamics (2.1) and that the ideal output y be identical 
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tc the model output y . Thus: 
r 

= Av (t) + Bf *(t) 

* * * * 

y (t) = Cv (t) ; v (o) = v 

o 

* 

where the ideal state v (t) is (as with v(t)) 
Hilbert space H. 


(3.1a) 

(3.1b) 


in the infinite dimensional 


Furthermore 

y*(t) - y m (t) = c m q(t) (3.2) 

In a manner similar to that in [5], it will be assumed that v (t) and 

* 

f (t) are linearly related to the model state vector q(t) and command vector 

u (t) as follows: 
m 


V (t) 
* 

' A 11 q(t) + S 12 u m 

(3.3) 

f (t) 

‘ S 21 q<t) + S 22 u m 

(3.4) 


The bounded linear operators S^, S 19 , will not be determined to 

satisfy MFC 1. 


To this effect, differentiation of (3.3) with respect to t and substitu- 
tion of (3.1) and (2.18) gives: 


9v (t) 
at 



11 


A q + S 
m 


11 


B u 
m m 


* 

Av 


+ Bf 


* 


(3.5a) 


where 


v = q + S. u 
o 11 n o 12 m 

is in D(A). 

Replacing v and f on the right side of (3.5) by (3.3) and (3.4) gives: 


(3.5b) 


s n V : £ n V. 

A < S ll1 + S 12 "»> + B < S 21 q + S 22 V (3 - 6) 

Now since (3.6) must be valid for all q and u m> it is necessary that: 

S.. A * AS,. + BS.. (3.7) 

11 m 11 21 

S 11 B » ’ AS 12 + BS 22 (3 - 8) 
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(3.9) 


Finally the incorporation of (3.2) yields 
y*(t) = CS n q + CS 12 U m = y m = C^q 

Thus : 


CS 11 = C m 

(3.10) 

CS 12 = 0 

(3.11) 


In summary then eqs. (3.7), (3.8), (3.10) and (3.11) must be^solved in 
order to find S 0 ^ and S ^ which in turn define the ideal control f of Eq. (3.4), 

Recall however, that both MFC 1 and MFC 2 must both be satisfied. In 
order to satisfy MFC 2, it is useful to consider the equation for the error 

e * v* - v (3.1?) 

* 

which is in D(A) when v and v are both in D(A). Differentiation of (3.12) 
with respect to time gives: 

* 

3e _ dv 
at " at 

= Av + Bf - (Av* + Bf*) 

= Ae + B(f - f*) (3.13) 

This equation suggests that the actual model following control f be defined as: 
f - f* + G(y - y m ) 

•k it 

= f + G C(v - v ) 

= f* + G C e 0.14) 

Substitution of (3.14) into (3.13) gives: 

e = (A + B G C)e (3.15) 

P P 

where G:R ->R is a bounded linear operator. Thus if G is chosen such that 
(A + B G C) generates an exponentially stable C -semigroup, then the control f 
as defined by (3.14) will satisfy the conditions for model following. 

It i9 important to note that this controller is clearly finite dimensional. 
For implementation it is only necessary to "build" a finite dimensional refer- 
ence model and form the proper linear combination of its state vector and command 
vector. The gain operator G is also finite dimensional and should be chosen 
such that the decay of any transient caused by initial plant model output error 
is sufficiently fast. We summarize the above discussion as 
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Theorem 1 : If (A,B,C) is exponentially output stabilizable and there exist 

bounded linear operators S ll* S 12 * $21* and S ^2 such that (3.7) - (3.8) and 

(3.10) - (3.11) are satisfied, then the model following control (3.4) an£ (3.14) 

satisfies the model |ollowirg conditions MFC (1) and (2) and lim [v(t)-v (t)] = 

0 when both v and v belong to D(A). t-* 30 

o o 

FroT" [10], ue see that (A,B,C) is exponentially output stabilizable if and 
only if 11^ = N(C) X and H R = N(C) form a pair of stabilizing subspaces for (A,B). 

Note that dim = P which is the number of sensors (or actuators) used. The 

conditions for existence of the ideal trajectories (3.1) will be developed in 
the next subsection. 

3.2 Existence of Ideal Trajectories 

* 

The existence of ideal trajectories v (t) for the DPS (2.1) is determined 

by solutions S. 4 to the operator equations (3.7) * (3.8) and (3.10) - (3.11). 

1 J 

These can be rewritten as 


A 

B 

s n 

S 12 

C 

0 

- S 21 

S 22- 



“ S ll 

0 


A 

m 

B 

m 


0 

I_ 


C 

m 

0 


(3.16) 


N P 

where S - :R -*D(A) and S ’.R -*D(A) are bounded operators with finite-rank and 
N 1 p P P iZ 

$ 2 i : ^ "* R and S ">2 :R "* R are matT ^ ces appropriate dimension. Note that 

(3.16) describes a kind of aggregation (in the sense of Aoki) for the infinite- 
dimensional system (2.1) injo a finite-dimensional system (2.17). The existence 
of the ideal trajectories v (t) in (3.1) guarantees that such an aggregation is 
possible, i.e. the DPS (2.1) generates the ideal trajectories which correspond 
to those of the finite-dimensional model (2.18). 


In most situations, the ideal initial condition will be v =0; hence, 

from (3.5b) we would choose q =0 and S 10 “ 0, which correctly corresponds to 

o Iz 

(3.11). This reduces the other operator equations to the following: 


S 11 A m A S n + B s 21 

(3.17a) 

S 11 B m " B S 22 

(3.17b) 

C S 11 = C m 

(3.17c) 


we have the following: 


Theorem 2 : If the spectra c(A) and a (A ) are separated by a smooth simple 

closed curve r containing o(A ) in its in?arior and a (A) in its exterior, then, 

m N P 

given any linear operator S-. :R -*-R , there exists a unique bounded linear 
N 

operator S^:R -*D(A) given by 
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(3.18) 


S ll q = 2^1 / R(X ’ A) B S 21 R(X »V qdX 


lor any q m R 


PROOF: From (3.17a), it follows that for any Xea(A)0 a (A ): 

m 

S 1t R(A,A ) - R(A,A) B S_ R(A,A ) - R(A,A) S 1t 
ii m 21 m 11 

But integration of (3.19) over the curve T produces: 


(3.19) 


0 - 2^1 J R ( X > A ) : n T dX " 27f / lS ll R < X »V q ' R(X * A) B S 21 R 0 *V qldX 


= S xl q - I R(A,A) B S„ R(A,A rt )q dA 


2-rri 


21 


because r encloses the finite number of singularities of A and excludes all of 
the spectrum of A. Clearly, since R(A,A) :H-*-D(A) , must°have its range in 

D(A), and this is the desired result. # 

Once, we have specified the matrix the unique operator is deter- 

mined. Satisfaction of (3.17c) could most easily be done by defining C to be 

m 

C S^. The determination of the matrix S ^2 for (3.17b) could be done from 

S._ = (B*B) _1 B* S., B (3.20) 

22 Ij. m 


as long as B is c sen so that a solution exists. Note that the operator B has 
10 * 
full ran* P and so the inverse of B B exists. 


Although the above existence result does not really require the number of 
actuators and sensors to be equal, this will be needed in the later sections. 
Also, the following alternative existence result requires it: 


-1 P 

Theorem 3 : Let zero belong to p (A) and C A B be nonsingular on R , then 


A B 

-l ] 

"®H n !2 _ 


A _1 (I-B(CA" 1 B) -1 CA“ 1 ) • A -1 B(CA“ 1 B)" i 

c o] 


_«2i ^ 22 - 


| 

_ (CA“ 1 B)" 1 CA“ 1 j -(CA _1 B) _1 


and S 


'll 


12 

satisfies: 


fi ll S 11 V S 21 


°21 S 11 A m + 


n 22 C m’ and S 22 


“21 S 11 ®m Where 


11 


n,, s.. a + ft. 9 c 

11 11 m 12 m 


(3.21) 


The proof of Theo. 3 can be obtained by straightforward computation using (3.16). 
Furthermore, note that 
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AS 


11 


Aft 


11 


S 11 A m + ^12 


C 

m 


- <i-«n 12 > 


S U A „ + <- Bn 22> C » 


■ S 11 A „ - B l°12 S U A * + °22 C J 


S 11 A m " BS ;i 


which is the same as (3.17a); however, Theo. 3 gives a wider range of solutions 
than Theo. 2 since S^ 2 need not be zero . The solution of (3.21) can be handled 

when zero belongs to (A ) because we then have the following: 


s n A 
11 m 


-1 


= ft 


11 11 


+ ft 


12 


C A 
m m 


-1 


(3.22) 


which has a unique solution whenever the o(A m ^) and a(fi.^) are separated 
by a smooth simple closed curve (see proof of Theo. 2). 


4. THE ADAPTIVE MODEL FOLLOWING CONTROLLER 


4.1 Development of the Adaptive Controller 

The nonadaptive control law (3.14) requires exact knowledge of the gain 

operators G, S ?1 , and S ?2 * These may be known to exist via mathematical structure 

of the DPS (A,B,C) in (271) (e.g. Theos. 1, 2, 3) but they may not be available 

in an explicit form. Consequently, we would need an adaptive version of (3.14): 

f(t) = S n (t) q (t) + S 22 (t) u m + G(t) e y (t) (4.1) 


where 



We assume throughout Sec. 4.0 thaj the hypotheses of Theo. 1 are satisfied 
for the DPS (2.1). Take e(t)* = v(t) - v (t) and, from (2.1), (3.1), (3.3) and 
(4.2), obtain (for v q and V q in D(A)): 


where 


Se(t) . 
8t 

e (o) 5 


A e(t') 4 BAK(t) 

c 


e 

o 


v 

o 


* 

V 

o 


r(t) 


A r A + BGC generates an exponentially stable C -semigroup U (t) and 


r(t) 


e y (t) 

q(t) 


u 


Ml Op 

belongs to R and AK(t) = K(t) - K q where 


m 
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K(t) - [ G ( t ) j S (t) j S 22 (t) ] and K - [G j S . J S ] 

I I * I 

The adaptive gain laws we sh*ll use ar: motivated by [6] ana have the form: 

K(t) * Kj(t) + K p u) (i 

K p (t)z - -r p e y (t) (r (t) ,z) (/ 

. K x (t)z = -Tj" 1 e y (t) (r(t),z) (t 

• _ N+2P 

where , z belongs to R , and T p , are both positive definite 

matrices on R . Note that (since K q is constant): 


AKj(t) = Kj(t) = -r x e (t) (r(t) , •) 


(4.5) 


where 


AK x (t) = r. z (t) - K q . 


The closed-loop adaptively controlled DPS is given by (4.3) and (4.5): 


9e(t) 


= A c e(t) + F (t, e(t)) 


e(o) = e = 
o 


K x (o) 


(4.6) 


where 


e(t) = 


F(t,e(t)) = 


AK x (t) 


BAK(t) r (t) 


A 0 
A = C 

c L° ° 


, and 


e y (t) 


, with e (t) * C e(t) and r(t) « o 

: r i «/t) (t(t),-)J >' [ u 


The state e(t) of (4.6) resides in a new Hilbert space H where H = H x B (k , 

P ^ 

R ) with B 2 (H if H 2 ) representing the Schmidt class of compact linear operators 

from into Hj with inner product (A,B) = tr A E where "tr" denotes the trace 

of the operator; see [11] pp 262-264 for further details. The inner product on 
H is formed by summing those of H and B 7 ; we shall use the same symbols for ail 
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inner_products_( • , •) and their corresponding norms j|*||. The nonlinear func- 
tion F(t,*)‘fr*H is continuous; hence, 

e(c) = U(t) e Q ; t ^ 0 (4.7) 

where U(t) is the nonlinear semigroup defined on 11 by (for any h in H) : 

U(t)h = U^(t)h + J U c (t-t) F(T„U(T)h)dr (4.8) 

where _ 

[U (t) 0 

U (t) = is the linear C -semigroup generated on H by A in 

c Lo ii 0 c 

(4*6), The above follows from [12] Lensna 5.2 p. 186 where further details on 
nonlinear semigroups are also available; consequently, the closed-loop infinite- 
dimensional system (4.6) xs well-posed on H. 

4.2 Closed-Loop Stability 

The stability analysis of the nonlinear infinite-dimensional system (4.6) 
requires the extension of Lyapunov theory to infinite-dimensional spaces. This 
has been done in [12] -[13] and we summarize the necessary elements here: 

Def : The equilibrium point $_is stable for the system _(4.6) if for every e > 0 

there exists t > 0 such that j je(o) - VJ \ < <5 implies ] j e (t) - t for all 

t ^ 0. If, in addition to stability, there is a y > 0 such that j |e(o) - <J> | | < 

Y implies lim | |e(t)-$| | = 0, then ♦ is said to be asymptotically stable for 

t'X® 

(4.6). Usually we can take 4=0. We say an equilibrium point is unstable 
whenever it is not stable. 

"?ef : t A continuous functional V:H-*R is a Lyapunov function for (4.6) if V(o) = 0 
and V(e) £ 0 for all e in H wnere 

V(e) = lim sup (4. 9) 

t-K) + 

where e is in H and e(t) = U(t)7 a 9 given in (4.7). 

Lemma 1 : If V:H+R is a Lyapunov function for (4.6) with the property that 

vG) 1 f 1 (i|e||) (4.10) 

for all e such that ||e|| ^h (where 0 < h < “) and f^ is of class (i.e. 
f ,:(0,taj->-R + with f^(o) = 0 and f^ strictly increasing on [0,h], then the zero 
equilibrium point is stable for (4.6). 

Lemrar 2 : If in addition to the hypotheses of Lemma 1, the Lyapunov function V(*) 

has th property: 
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V(e) <_ -W(e) for all e in H 
W(e) >_ f ( 1 ! e j j ) for ||e|| <^h 


(4.11b) 


where is also of class M. , then the zero equilibriua point is asymptotically 
stable hr (4.6). * 

The proofs of Lemmae 1 and 2 can be found in fl3J. These results constitute 
Lyapunov's Direct Method on infinite-dinenaional spaces. 

We now have the following stability result for our adaptively controlled 
closed-loop system (4.6): 

Theorem 4 : Assume the following: 


(a) In (4.3), A c = A + BGC satisfies 

(A v, Pv) + (Pv, A v) * -(Qv,v) 

c c 


(4.12) 


for all v in D(A) where P and Q are symmetric positive operators on H such that 
(for some ci, S positive constants): 


I I v | | 2 1 (v,Pv) <6 | |v| | Z 

u|lv| | 2 <_ (Qv,v) (i.e Q is coercive) 


(4.13a) 

(4.13b) 


for all v in H, 
(b) B*P = C, 


(4.14) 


(c) the hypotheses for Theo. 1 are satisfied, and both v q and v q belong to D(A) , 


then V(e) = (e,Pe) + (AKj, Tj AK^)» with AK^t) 5 K^t) - K q and e = |^|» is a 
Lyapunov function for (4.6) and the zero equilibrium point is stable. 


PROOF : 


Recall that 


AK(t) = AK^t) + K p (t) 
AK I (t) = K x (t) 


(4.15) 

(4.16) 


Now, clearly V is a continuous functional from H into R (due to (4.13a) with 
V(o) * 0. Furthermore, since V is a quadratic functional, it is Frechet 
differentiable. Hence, from (4.6) and (4.12), 

V(e) = -(Qe,e) + 2p (4.17) 

where p = [(Pe, BAKr) + (AK^.F^ A*K^)] 

From (4.16), (4.4c), and (4.15), we have 
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p = (B*Pe, AK r ) - (AKj, e y (r,-)) 

= (B*Pe , AK r ) - (r, AK^ e y ) 

= <B*Pe, A^r) + (B*Pe, K p r) - (r, AK* e y ) 

= (AKjr, [B*Pe - e y ]) + (K p r, B*Pe) (4.18) 

* * 

where we have used (A,B) = tr A B = tr(BA ). Furthermore, using (4.14) in 
(4.18), yields 

V « (K p r, e y ) = -(r p e y )||r || 2 (4.19) 

from (4.4b). Consequently, using (4.19) in (4.17), we obtain 
»W‘ -KQe.e) + 2(r p e y . e y )l|rl| 2 l 

I- MM I 2 ♦ 2 « p l!e y |! 2 ||r|| 2 ) < 0 (*.20 

where a = X , (T ) and we have used (4.13b). 
p min p 

Also, using (4.13a), we have 
V(e) > | | e | | 2 + A nin (r i ) | | AK X | j 2 

In other words, f, ( 4 ) = [1 + X . (r 
l min 

abo^e satisfies the hypotheses of Lemma 1 and the desired result is obtained. # 

Note that the use of a proportional adaptive gain (4.4b) produced the second 
term in (4.20): however, this term is not essential and the above argument 
could be simplified by omitting (4.4b) from the adaptive gain laws. 

The hypotheses (a) and (b) correspond to the Kalman - Yakubovich conditions 
in infinite-dimensional spaces. From [13] Theo. 4.7, if for some real w, 

(Av, v) u>| | v | | for all v in D(A), then exponential output stabilization of 

(A,B»C)would be equivalent to satisfaction of hypotheses (2); however, there 
would be no guarantee that ? and 0 could be found in (4.12). such that (4.14) 
could be obtained. In finite-dimensional spaces, the Kalman - Yakubovich 
conditions are equivalent to the strict positive realness of the transfer 

function T c (s) = C(sl-A c ) ^B, i.e. Re T c (ju) > 0 for all real 01 ; see [14] pp. 

115-118. A number of papers, e.g. [15] - [17], have been written on this 

relationship in infinite-dimensional spaces. For example, [17] asserts that 

ReT (j<*>) must be coercive*which would be quite a bit stronger than what is 
c 

required in finite-dimensions. This is an area that requires further investi- 
gation. 

As pointed out in [9], we cannot immediately conclude asymptotic stability 
from (4.20) since it does not satisfy the hypotheses of Lemma 2. In finite- 
dimensional space, we could apply the LaSalle Invariance Principle to obtain 


which is of class 


Therefore, the 
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asymptotic stability as is done in [6]; however, in infinite-dimensional spaces, 
it is not the case that '’bounded sets are precoapact** and this is essential for 
the LaSalle result. 


The following result ([13] Theo. 5.4 p. 188) may be helpful: 

L emma 3 : Let in (4.6) generate the linear C^-semigroup U^(t) on E and F 

is any_bounded x continuous function such that (4.6) generates a nonlinear semi- 
group U(t) on H (as given in (4.8)), then all bounded orbits of (4.6) are pre- 
compact if either 


(a) U q ( t) is compact operator for all t 0 
or 


(b) U c (t) is exponentially stable and the function F is compact (i.e. maps 

bounded sets into precoopact ones) 

A 0 

— _ c 

Due to the form of A = , it is not possible to satisfy (b); however, 

c [0 0 

(a) may be satisfied, for example by operators A which generate holomorphic 
semigroups. This latter is determined by the form of damping operator in a 
flexible structure. Again, this is a topic for further investigation. An 
alternative adaptive gain lav: 


KjfOv = -r 


-1 


(e y (r,v) + K I (t)v) 


(4.21) 


yields : 

V(e)< -[a | | e | | 2 + 21'iAKjjl 2 + 2( AKj, K q )] 


which does not quite give asymptotic stability but might be modified to do so. 

5 . CONCLUSIONS 

In this paper, we have presented a direct adaptive controller for linear 
distributed parameter systems (DPS) described on infinite-dimensional Hilbert 
spaces. The controller is based on a command generator tracker approach used 
in finite-dimensional spaces, e.g. [6] where it is shorn to be asymptotically 
stable. We have shown here that, under certain conditions on the open-loop 
loop DPS, ideal trajectories do exist and the adaptive controller is stable, 
i.e. the output and gain errors remain bounded. If the further condition that 
A in (2.1) generates a holomorphic C -semigroup is imposed, then we can also 
conclude asymptotic stability which guarantees asymptotic tracking or model 
following. 

A number of issues have been opened for further investigation: 

(1) use of dynamic rather than output feedback stabilization; 
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(2) generation of asymptotic ideal trajectories by the open-loop DPS; 

(3) connections between the Kalman-Yakubovich conditions and the input- 
output description of the DPS; 

(4) development of alternative adaptive gain lavs which produce asymptotic 
stability of the closed-loop system; 

(5) exploration of reasonable conditions under which LaSalle's Invariance 
Principle can be used to determine asymptotic stability of the closed- 
loop system. 
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